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ABSTRACT 

Recent observations have shown that in many exoplanetary systems the spin axis of 
the parent star is misaligned with the planet's orbital axis. These have been used 
to argue against the scenario that short-period planets migrated to their present-day 
locations due to tidal interactions with their natal discs. However, this interpretation is 
based on the assumption that the spins of young stars are parallel to the rotation axes 
of protostellar discs around them. We show that the interaction between a magnetic 
star and its circumstellar disc can (although not always) have the effect of pushing the 
stellar spin axis away from the disc angular momentum axis toward the perpendicular 
state and even the retrograde state. Planets formed in the disc may therefore have their 
orbital axes misaligned with the stellar spin axis, even before any additional planet- 
planet scatterings or Kozai interactions take place. In general, magnetosphere-disc 
interactions lead to a broad distribution of the spin-orbit angles, with some systems 
aligned and other systems misaligned. 

Key words: accretion, accretion discs - planetary systems: protoplanetary discs - 
stars: magnetic fields 



1 INTRODUCTION 
1.1 Background 

Transiting planets are providing new ways to charac- 
terize exoplanetary systems. In particular, the Rossiter- 
McLaughlin (RM) effect, an apparent radial velocity 
anomaly caused by the partial eclipse of a rotating parent 
star by its transiting planet, can be used to measure the 
sky-projected stellar obliquity, the angle between the stellar 
spin axis and the planetary orbital axis. As of August 2010, 
sky-projected stellar obliquity has been measured in 28 sys- 
tems fusing the RM effect (see Triaud et al. 2010; Winn 
et al. 2010). Among these, about 60% have an orbital axis 
aligned (in sky projection) with the stellar spin, while the 
other systems show a significant spin-orbit misalignment, 
including 5 systems with retrograde orbits (e.g., Hebrard 
et al. 2008: Winn et al. 2009; Johnson et al. 2009; Narita 
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1 The number has increased to 48 in October 2010 (A. Triaud, 
private communication). 



et al. 2009; Pont et al. 2010; Jenkins et al. 2010; Triaud et 
al. 2010). In addition, a recent analysis of the stellar rotation 
velocity shows that in 10 out of a sample of 75 exoplanetary 
systems there is likely a significant degree of misalignment 
along the line of sight between the planetary orbital axis and 
the stellar spin axis (Schlaufman 2010). 

It is generally accepted that close-in exoplanets ("hot 
Jupiters") are formed at a distance of order several AU's or 
larger from their host stars before migrating inwards to their 
current locations. Gravitational tidal interaction between a 
gaseous disc and a young planet (Goldreich & Tremainc 
1980) provides a natural mechanism for the inward planet 
migration (Lin et al. 1996; see Lubow & Ida 2010 for a re- 
cent review). Gas-driven migration alone, however, is un- 
likely to explain the observed eccentricity distribution for 
planets with periods longer than a few weeks. Two other 
processes have been suggested to play a role in shaping 
the architecture of exoplanetary systems: (i) Strong gravita- 
tional scatterings between planets in a multiplanet system 
undergoing dynamical instability (e.g., Rasio & Ford 1996; 
Weidenschilling & Marzari 1996; Zhou et al. 2007; Chat- 
terjee et al. 2008); (ii) Secular Kozai interactions between a 
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planet and a distant companion (a star or planet) in a highly 
inclined orbit (Wu & Murray 2003; Fabrycky & Tremaine 
2007; Wu et al. 2007; see also Eggleton & Kiseleva-Eggleton 
2001) or between planets in a multiple systems (Nagasawa 
et al. 2008). Both scatterings and Kozai interactions are ex- 
pected to cause eccentricity in the final planetary orbits and 
misalignment between the stellar spin and the planetary or- 
bit axis. 

Newly formed planets are expected to lie in same plane 
as the gas disc. It is usually presumed that the planet-disc 
interaction preserves the alignment between the gaseous disc 
and the orbit of the planet. Under the assumption that the 
stellar spin is aligned with the disc angular momentum axis, 
planetary systems with zero or small stellar obliquity would 
be produced. The discovery of a significant fraction of mis- 
aligned systems (particularly the retrograde systems) has 
been suggested as a blow against the theory of disc-driven 
migration, and to favor the Kozai cycles plus tidal interac- 
tions as the primary mechanism for the formation of hot 
Jupiters (e.g., Triaud et al. 2010; Winn et al. 2010). 

Planets formed at a few AU's are scattered to the prox- 
imity of their host stars when their eccentricities approach 
unity. Their orbits may be circularized by subsequent tidal 
dissipation in the planets but its efficiency depends sensi- 
tively on the distance of closest approach (e.g., Ivanov & 
Papaloizou 2004). Dynamical relaxation generally leads to 
a Rayleigh distribution in eccentricity (Zhou et al. 2007). If 
this process is the leading cause for planets to venture into 
their stellar proximity, its combined effect with the tidal 
circularization process would yield a continuous periastron 
distribution (Zhang et al. 2010, in preparation), which is not 
consistent with the observed sharply bimodal distribution. 

Current data suggest that there may be two popula- 
tions of short-period exoplanet systems, one with spin-orbit 
alignment, the other with significant misalignment (Schlauf- 
man 2010; Winn et al. 2010). Although stellar or planetary 
companions have been identified as potential culprits for 
Kozai mechanism to operate in some systems, it may not be 
the dominant process to account for the origin of many hot 
Jupiters which do not appear to be either bound to binary 
stars or associated with planetary siblings with comparable 
masses and periods less than a decade. 

The solar system also provides a clue. Except for Pluto, 
all planets outside 1 AU lie within 2° of the ecliptic plane, 
while the Sun's equatorial plane is inclined by 7° with re- 
spect to the ecliptic. There is no obvious celestial candidate 
which can impose sufficient secular perturbation to induce 
this observed spin-orbit inclination. 

We note that spin-orbit misalignment is not limited to 
planetary systems. Hale (1994) has measured the inclination 
to the line of sight of the spins of stars in binaries by compar- 
ing the rotational period to the v sin i values of rotationally 
broadened lines, and inferred that binaries are spin-aligned 
for a <y 30 — 40 AU, but become randomly oriented for larger 
orbits. A particularly striking example is the binary system 
DI Herculis: Both stars rotate with their spin axes nearly 
perpendicular to the orbital axis and inclined to each other 
(Albrecht et al. 2009). Recently, a number of close binaries 
(with period of a few days) have been found to have nonzero 
stellar obliquities (A. Triaud, J. Winn, private communica- 
tions, 2010). 




Figure 1. A sketch of the effect of the magnetic warping torque. 
This torque tends to push the disc angular momentum axis I 
toward a state perpendicular to the spin axis ui a . When I is fixed 
by the outer disc, the back reaction torque tends to push the spin 
axis toward being misaligned with I. 

1.2 This Paper 

In this paper and the companion paper (Foucart & Lai 2010, 
hereafter Paper II), we explore an alternative scenario for 
the observed spin-orbit misalignment in exoplanet and bi- 
nary star systems. We assume that the planets' present-day 
orbits are in the plane of their natal discs. We study the 
mutual interaction between a magnetic young star and its 
surrounding disc, and show that under certain (but real- 
istic) conditions, the stellar spin axis can be pushed away 
from the disc axis toward the perpendicular state and even 
the retrograde state. Planets formed in the disc therefore 
may have their orbital axes misaligned with the stelar spin 
axis, even before any additional planet-planet scatterings or 
Kozai interactions take place. The basic idea is the follow- 
ing. A magnetic protostar (with B* <; 10 3 G, typical for 
classical T Tauri stars; see, e.g., Bouvier et al. 2007; Donati 
& Landstreet 2009) generally exerts a warping torque and a 
precessional torque on the inner region disc before the disc 
is disrupted at the magnetosphere boundary (Section 2; see 
Lai 1999). These torques have a tendency to make the inner 
disc tilt away from the stellar spin axis and precess around 
it on a warping timescale, which is much longer than the 
rotation period (see Fig. [3] below for a simple "laboratory" 
toy model that explains the origin of the warping torque) 
Q However, internal processes in the disc will try to resist 
the inner disc warping, either by viscous stress or by bend- 
ing wave propagation. The result is that, for the reasonable 
disc/stellar parameters, the inner disc may not be signifi- 
cantly warped (see Paper II), i.e., the direction of the inner 
disc is approximately aligned with that of the outer disc (as- 
sumed to be fixed; but see below). The back-reaction torque 

2 Note that the secular disc warp discussed in this paper is dif- 
ferent from the dynamical warp which varies on the timescale of 
the stellar rotation period. Such dynamical warp may arise from 
the periodic vertical force on the disc from the magnetic star 
(Terquem & Papaloizou 2000; Lai & Zhang 2008), or from the 
simple effect where the disc material "climbs" up the field lines 
at the magnetosphere boundary before funneling to the magnetic 
polar cap [see Bouvier et al. (2007) for possible observational evi- 
dence of dynamical warps in the classical T Tauri star AA Tauri] . 
The dynamical disc warp averages to zero over an rotation period 
and has no effect on the secular evolution of the system. 
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will then act on the star, changing the spin direction on a 
timescale much longer than the disc warping time. Thus, 
even if initially the stellar spin axis is approximately (but 
not perfectly) aligned with the disc axis, given enough time, 
the stellar spin axis will evolve towards the perpendicular 
state and even the retrograde state (see Fig. [TJ. Therefore, 
a young planet, formed with its orbit axis aligned with disc 
axis, may be misaligned with the stellar spin axis. 

While the magnetic warping torque alone pushes the 
stellar spin towards misalignment, other torques (such as 
that due to the angular momentum carried by the accret- 
ing gas) tend to align the stellar spin with the disc axis. 
When these torques are included, the misalignment effect 
is reduced. However, we show that with reasonable disc pa- 
rameters, the warping torque could still be dominant and 
spin-disc misalignment can be (but not always) produced 
under general conditions. 

We note that even without the magnetic effect discussed 
above, misalignment between stellar spin and disc axis could 
be a natural consequence of star formation process itself. 
Indeed, when stars form under the turbulent conditions of 
molecular clouds, the gas accreting onto the protostar does 
not necessarily fall in with a fixed orientation (e.g., Mc- 
Kee & Ostriker 2007), i.e., the outer disc direction can vary 
in time. Additionally, the disc can be perturbed by other 
nearby stars and change orientation (e.g., Bate et al. 2003; 
Pfalzncr ct al. 2005). If the stellar spin is determined by 
the total angular momentum gained over its whole forma- 
tion history, while the orbit of a planet is coplanar with 
the accretion disc towards the end of its evolution, then 
spin-orbit misaligned planetary systems may be produced 
(Bate et al. 2010). However, since the gas disc carries a large 
amount of angular momentum, one might expect that effi- 
cient spin-orbit alignment will be achieved in the absence of 
any magnetic interaction. Thus, in this picture, it is impor- 
tant to understand how the stellar spin evolves and on what 
timescale (and whether this timescale is shorter or longer 
than the timescale of varying outer disc orientation). 

The remainder of this paper is organized as follows. In 
Sect. [2] we present an analytical model used to describe the 
interaction between a magnetic star and its disc, and derive 
the magnetic torques acting on the disc. This is similar to 
the model considered in Lai (1999), but takes into account of 
more recent works on magnetosphere - disc interactions. In 
Sect. 3 we qualitatively discuss disc warp due to the mag- 
netic torques, relegating technical details to Paper II. We 
then study in Sect. 4 the evolution of stellar spin axis due to 
the back-reactions of the magnetic warping torque and other 
torques, under the assumption that the disc is approximately 
flat. In Sect. 5 we apply our theory to the problem of spin- 
orbit misalignment in exoplanetary systems and present two 
possible scenarios that may play a role in explaining the ob- 
servations. Finally we discuss the implications of our results 
in Sect. 6. 



2 ANALYTIC MODEL OF MAGNETIC STAR - 
DISC INTERACTION: DISC WARPING 
TORQUE 

The interaction between a magnetic star and its circumstel- 
lar disc has been studied over many decades, both theoret- 



ically (e.g., Pringle & Rees 1972; Ghosh & Lamb 1979; Aly 
1980; Wang 1987; Aly & Kuijpers 1990; Shu et al. 1994; van 
Ballegooijen 1994; Lovelace et al. 1995,1999; Lai 1999; Uz- 
densky et al. 2002; Pfeiffer & Lai 2004; D'Angelo & Spruit 
2010) and using numerical simulations (e.g., Hayashi et 
al. 1996,2000; Miller & Stone 1997; Goodson et al. 1997; 
Fendt & Elstner 2000; Matt et al. 2002; Romanova et 
al. 2003,2009). Most previous works deal with the case where 
the stellar spin axis, the magnetic axis and the disc axis are 
aligned. Even for such a "simple" case, the problem is com- 
plex. There are at least four physical ingredients involved in 
the magnetosphere - disc interaction: (i) The stellar mag- 
netic field penetrates through part of the disc, establishing 
magnetic linkage between the star and the disc. This may 
be achieved either through dissipation in the disc (if the disc 
is sufficiently dissipative) or, more likely, through reconnec- 
tion between the stellar and disc fields - the latter may be 
associated with dynamo actions in the disc, (ii) In the in- 
ner disc regions where ionization fraction is non-negligible, 
magnetic fields diffuse through the disk on a time scale much 
longer than the orbital period (e.g., Shu et al. 2008; Terquem 
2008). In this case, the field lines linking the star and the 
disc are twisted because of the difference in the stellar ro- 
tation oj s and disc rotation f2(r), generating toroidal fields 
Ai?0 from the vertical field B z that threads the disc. Af- 
ter a characteristic time of order \uj b — f2(r)| _1 , the toroidal 
field |AB^| becomes comparable to \B Z \, and the flux tube 
starts expanding at a fast rate and the fields open up, tem- 
porarily disconnecting the star - disc linkage (e.g., Aly 1985; 
Aly & Kuijpers 1990; van Ballegooijen 1994; Lynden-Bell & 
Boily 1994; Lovelace et al. 1995). However, reconnection be- 
tween the inflated field lines relaxes the shear and restore 
the linkage between the star and the disc, and the whole 
cycle repeats (Aly & Kuijpers 1990; van Ballegooijen 1994; 
Goodson et al. 1997; Matt et al. 2002^ (iii) Some open field 
lines are strongly "pinched" and twisted by the conducting 
disc, leading to a conical wind from a localized region near 
the magnetosphere-disc boundary which carries away angu- 
lar momentum from the accreting gas (e.g., Shu et al. 1994; 
Romanova et al. 2009). (iv) In the open field line region 
near the stellar rotation axis, outflows could be launched, 
carrying away angular momentum from the star (e.g., Matt 
& Pudritz 2005). Even without significant mass loss, angu- 
lar momentum may be transferred outwards through Alfven 
waves. 

All these four processes are likely to play a role in the 
spin evolution of the protostar. Obviously, for the general 
case where the stellar spin axis, magnetic axis and disc 
axis are misaligned, the situation is even more complicated. 
Moreover, the magnetic warping and spin-orbit misalign- 
ment considered in this paper take place on timescales much 
longer than the dynamical time (i.e., the spin period or disc 
orbital period), and therefore cannot be easily captured in 
3D numerical simulations. 

Nevertheless, the key physical effects of the magnetic 
star - disc interaction relevant to this paper can be de- 
scribed robustly in a parametrized manner (Lai 1999; see 



3 Strictly steady-state models have also been considered, but they 
generally require special conditions for the disc (see Uzdensky 
2004 and references therein). 
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Figure 2. A sketch of magnetic field configuration in a star — 
disc system for nonzero (the angle between the disc axis and 
the stellar spin axis) and 6+ (the angle between the stellar dipolc 
axis and the spin axis). Part of the stellar magnetic fields (dashed 
lines) penetrate the disc in the interaction zone between the disc 
inner radius n a and r lnt in a cyclic manner, while other field lines 
are screened out of the disc. The closed field lines are twisted by 
the differential rotation between the star and the disc, which leads 
to a magnetic braking torque and a warping torque. The screening 
current in the disc leads to a precessional torque. 



Fig. [2}. The stellar magnetic field disrupts the accretion flow 
at the magnetospheric boundary. Some of the accreting gas 
are channelled onto the polar caps of the star while other gas 
could be ejected in an outflow. The magnetosphere boundary 
is located where the magnetic and plasma stresses balance. 
For a dipole field, the inner disc radius is 
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where K'h and [i are the mass and magnetic dipole moment 
of the central star, M is the mass accretion rate, and 77 is a 
dimensionless constant of order unity, with typical estimates 
ranging from 0.5 to 1 and recent numerical simulation giving 
r\ ~ 0.5 (e.g., Long et al. 2005). Before being disrupted, 
the disc generally experiences magnetic torques from the 
star. The vertical (perpendicular to the disc) magnetic field 
produced by the stellar dipole is given by 
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where 9+ is the angle between the magnetic dipole axis ft 
and the spin axis u a , and fi is the angle between 61 s and 
the disc axis I (see Fig. [2)1 (note that for a warped disc, I 
depends on r). We assume that the static field component, 
= — (fi/r 3 ) cos 8i, cos ft, penetrates the disc in an "in- 
teraction zone", between r- m and r mt . As discussed before, 
this field is twisted by the differential rotation between the 
star and the disc, undergoing cycles of field inflation and re- 
connection. The differential rotation generates toroidal field, 
whose value at near the disc increases in time until it be- 
comes comparable to |-Bi s '|, at which point reconnection re- 
duces the field twist, and then the cycle repeats. We param- 
eterize the averaged value of the toroidal field generated by 
the twist by 



where the upper/lower sign refers to the value above/below 
the disc plane, and \C\ ~ 1 (e.g., Aly 1985; van Ballegooijen 
1994; Lynden-Bell & Boily 1994; Lovelace et al. 1995). The 
sign of £ is such that £ > for S7(r) — ui a cos /3 > 0, and £ < 
for fi(r) — uj a cos /3 < 0. Note that near the corotation radius, 
SI = uj s cos ft, the toroidal field may be limited by field diffu- 
sion inside the disc. Thus, \AB^/ B z \ ~ |f2 — uj s cos f3\r <5 1 
near the corotation radius, where r is the field dissipation 
timescale for the toroidal field in the disc. However, such 
region is very small and can be neglected, and we will adopt 
Eq. throughout this paper. The twisted toroidal field |3| 
implies a radial surface current K r — {c/2n)C,Bi s \ The ver- 
tical field Bz acts on K r , giving rise to an azimuthal force 
on the disc material and a (well-known) magnetic braking 
torque (per unit area) 
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2tt 
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This torque exists even for the aligned case. For j3 7^ 0, 
however, there exists an additional torque: The interaction 
between K r and the toroidal component of the dipole field, 
Si = — (/i/r 3 )(/x ■ cj>) (where <j> is the unit vector is the 
azimuthal direction), gives rise to a vertical force on the 
disc: 



F* = i [(B<"> + CM S) ) 2 - (B<? - (Bi 



2tt * z ' 
(5) 

After averaging over the azimuthal angle in the disc and 
the stellar rotation period, the net torque (per unit area) on 
the disc is 



N„ 



47T7" 5 



- COS /3 COS di, I X (u> 3 X I ) . 



(6) 



For £ > (i.e., inside the corotation radius), the effect of 
this torque (for a fixed spin axis o> s ) is to push the disc axis I 
away from £j s toward the "perpendicular" state (see Fig.[T]). 
The characteristic warping rate is (for (3 ~ 0) 



r„(r) 



Cm 2 



47rr 7 n(r)E(r) 



2 a 
COS 0* 
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where E(r) is the surface density of the disc. Note that in 
the case of /3 = 0* = 0, the corotation radius r co (where 
SI = uj s ) is somewhat larger than ri n for stars in the spin 
equilibrium (e.g., the simulation by Long et al. 2005 indi- 
cates that r co /ri n lies in the range of 1.2-1.5). Thus, for the 
inner region of the disc most relevant to our paper, the con- 
dition £ > is satisfied. 

The warping torque also exists for more complex stellar 
magnetic fields, which may be present for accreting T Tauri 
stars (e.g., Hussain et al. 2009; Donati et al. 2010). Consider 
an axisymmetric quadrupole field with the symmetric axis 
along Q (the 3rd axis) such that the magnetic quadrupole 
moment tensor satisfies Qn — Q22 = — Q33/2 = —Q/2. 
Outside the star, the stellar quadrupole field is given by 



3Q 
4r 4 



[5(Q-f) 2 -l] 



2r 4 ^ 



f)Q. 



(8) 



ABd 



(3) 



Again, we assume that the static component of the verti- 
cal field, B { q] z = - (3Q/4r 4 ) (3 cos 2 <9q - 1) sin cos /3 sin <j> 

(where 9q is the angle between a> s and Q), penetrates 
through the disc in a finite region between rj n and r- m %, and 
gets twisted to produce A_B^ = ^QBq z . After azimuthal 
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Figure 3. A toy model for understanding the origin of the warp- 
ing torque. A tilted rotating metal plate (with angular momentum 
L) in an external magnetic field B experiences a vertical magnetic 
force around region 2 and 4 due to the interaction between the 
induced current J and the external Bm, resulting in a torque N 
which further increases the tilt angle f). 



averaging and time averaging over the stellar rotation, we 
find the warping torque due to the quadrupole field: 



N 



(Q) 



1287rr' 



(3 cos 



-l) 2 sin 2 /3cos/3 Zx(w 3 xZ). (9) 



Thus, this torque has the same qualitative effect as the warp- 
ing torque due to the magnetic dipole. For simplicity, we will 
nee lect NL Q) in the remainder of this paper. 

We note that the existence of the warping torque is 
rather robust. The two basic ingredients are: (i) Magnetic 
field lines from the star penetrate the inner disc which ro- 
tates faster than the star; (ii) the connected field lines are 
twisted, reaching a quasi-steady twist angle. Both ingredi- 
ents involve dissipations. Qualitatively speaking, the warp- 
ing torque tends to push the system toward the perpendic- 
ular state (with the disc axis perpendicular to stellar spin) 
in order to minimize the energy associated with the twisted 
magnetic fields. In fact, the essential physics of the warping 
torque is captured by the following "laboratory" toy model 
(cf. Lai 2003), depicted in Fig. [3] Consider a metal plate 
(not a perfect conductor) in an external uniform magnetic 
field B, with the plate surface initially perpendicular to B. 
Neglect gravity. When the plate rotates (due to some exter- 
nal torque) with the rotation axis L, an EMF is induced (in 
the direction of v x B) and a radial current J generated in 
the plate - this is the same surface current K r mentioned 
before eq. Q. The interaction between J and B gives rise to 
a "magnetic braking" torque which tends to slow down the 
plate's rotation. Now tilt the plate by an angle /3 (and keep 
the plate rotating). The interaction between J and By (the 
component of B parallel to the plate) produces a vertical 
force around region 2 and region 4 in the plate. The result 
is a torque N which tends to push the plate's rotation axis 
L away from the magnetic direction, further increasing the 
tilt. 

In addition to the warping torque, there is also a pre- 
cessional torque on the disc when /3 / 0. This arises from 
the dielectric property of the disc: If the disc does not allow 
the vertical stellar field (e.g., the rapidly varying component 
of B z due to stellar rotation) to penetrate, an azimuthal 
screening current K$ will be induced in the disc. This K^, 



interacts with the radial magnetic field B r from the stellar 
dipole and produces a vertical force. After azimuthal aver- 
aging and averaging over the stellar rotation, we obtain the 
torque per unit area: Q 

,2 



N„ 



sin 8+ cos j3Cj s x Z, 



Tr 2 r 5 D(r) 

where the dimensionless function D(r) is given by 



(10) 



D(r) = max (jr^/r^ - 1, ^2H(r)/r in ^ , (11) 



and H (r) is the half-thickness of the disc. The torque ( 10 1 
tends to make the disc precess around the stellar spin axis. 



3 WARPED DISCS 

As shown in Sect. 2, the inner region of the disc where mag- 
netic field lines connect star and the disc and where the disc 
rotates faster than the star (SI > fi s cos /3) experiences a 
warping torque and a precessional torque. If we imagine di- 
viding the disc into many rings, and if each ring were allowed 
to behave independent of each other, it would be driven to- 
ward a perpendicular state (! 1 ii s ) and precess around the 
spin axis of the central star. Obviously, real protoplanetary 
discs do not behave as a collection of non-interacting rings: 
Hydrodynamic force and viscous force provide strong cou- 
plings between different rings. Depending on the physical 
condition of the disc, the dynamics of a warped disc may be 
driven by viscous diffusion (when the viscous a-parameter 
is greater than the dimensional disc thickness H/ r) or prop- 
agation of bending waves (when a £ H/r) (Papaloizou & 
Pringle 1983; Paper II). 

Thus, for a given outer disc direction Z out and stellar 
spin axis a> s , the disc will generally evolve into a warped 
state, with Z dependent on r. In particular, the orientation 
of the disc at the inner radius, Zi n = Z(n n ), generally differs 
from Z ut- For typical protostar parameters, M* ~ 1 Mq, 
R+ ~ 2R Q , fj, = B+R'i (with the surface field ~ 10 3 G), 
and mass accretion rate M ~ 1O _8 M0 yr _1 , the inner disc 
is located at a few stellar radii. From equation Q, we find 
that the timescale for the warp evolution is of order 



Mi, 



1/2 



8R P 



11/2 



R* J \lM e 
v \ „ ,-i 



10 g cm - 



(C cos ' 



(12) 



As shown in Paper II, under the combined actions of warp- 
ing/precessional torques and disc viscosity or bending waves, 
the disc will settle down into a steady state on a timescale 
that depends on disc viscosity and sound speed [see related 
work by Papaloizou & Terquem (1995) for discs in binary 
stars with inclined orbits.] In general, this timescale can be 
up to several orders of magnitude of T^ 1 evaluated at ri n . 
Nevertheless, it is much shorter than the disc lifetime and 
the timescale of the secular evolution of the stellar spin 



(see Eq. 14 below). Moreover, we show in Paper II that 



4 This assumes that only the spin-variable vertical field is 
screened out by the disc. If the vertical field is entirely screened 



out, then (—sin 0*) should be replaced by (2 cos 



■ sin 
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for most reasonable stellar/disc parameters, the steady-state 
disc warp is rather small because of efficient viscous damp- 
ing or propagation of bending waves. Thus, we will adopt 
the approximation l(r) ~ Z out in the remainder of this pa- 
per, leaving more detailed study of the effect of warped discs 
to Paper II. 



4 SECULAR EVOLUTION OF STELLAR SPIN 

We now consider the evolution of the stellar spin direction 
due to the back-reaction of the magnetic warping torque and 
other torques. Note that we are interested in the long-term 
evolution of the stellar spin. The specific angular momentum 
of the accreting gas at r- m is given by the Keplerian value 
\/GM+r 'in, so the characteristic accretion torque on the star 
is 

No = M(GM*n n ) 1/2 . (13) 
Thus the fiducial timescale for the spin evolution is 



4R, 



(14) 



where we have assumed the stellar spin angular momentum 
J s = 0.2M*i? 2 cj s (for a fully convective protostar modeled 
as a r = 5/3 polytrope). 

In general, the evolution equation for the spin angular 
momentum of the star, J s a> s , can be written in the form 



d 
tit 



(J s u> s ) = Af = Afi+Af s + Af w + Af P 



(15) 



Here JV ; represents the torque component that is aligned 
with the inner disc axis l ln (= I for a flat disc) : 



Afi = AAf(GAf*r in ) 1/2 Z 



(16) 



where A is a dimensionless parameter. Equation ( 16 1 in- 



cludes not only the accretion torque carried by the accreting 
gas onto the star, M aec (GMriii) 1 ~ I (note that in general, 
the accretion rate onto the star, M acc , may be smaller than 
M, the disc accretion rate), but also includes the magnetic 
braking torque associated with the disc - star linkage, as 
well as any angular momentum carried away by the wind 
from the magnetosphere boundary. While the details are 
complex, we would expect that all these contributions to 
Af j tend to make A < 1, as suggested by recent numerical 
simulations (e.g., Romanova et al. 2009) and earlier works 
(e.g., Shu et al. 1994). The term Af s = — [jV s |w s represents 
a spindown torque carried by a wind /jet from the open field 
line region of the star. In the aligned case (I = a> s ) stud- 
ied by previous works, Afi and Afs are the only torques 
acting on the star, and spin equilibrium is reached when 
AM(GM*r in ) 1/2 = |JV S |. Note that in general, the values of 
A and N s may depend on P and other quantities. 

The term Afw and Af p represent the back-reactions 
of the warping and precessional torques, respectively. From 



Eq. (§, we have 



27rrN m dr 



a' 2 
C A* 



cos 9* cos pi x (lJ s x t) 



6rf n 

= Non w I x (w s x I) 



where 



witfJ3 



n w = /2 cos 2 6»* cos/3, 



C = C [1 - (n n /r int ) 3 ] 



Similarly, from Eq. (10 1, we have 



Af P = 
with (for thin discs) |^j 



: NoTlp lJ b X i, 



cos p. 



(17) 



(18) 



(19) 



(20) 



(21) 



Note that both Af w and Afp are of order /x 2 /rf n (see the 



first line of Eq. 



17 



which does not directly depend on M. 



However, when we use Eq. for Hn, we find that both Af u 
and Afp are of the same order of magnitude as the fiducial 
accretion torque No. Therefore, as we will see below, the 
timescale to change the stellar spin direction (if at all) is of 
the same order as t SD j n . 



From Eq. ( 15 1, we find that the magnitude of the stellar 



spin evolves according to 



d_ 
dt 



J s =Af- 



— No (A cos p + n w sin 2 /3) — |7V S 



(22) 



The inclination angle of the stellar spin relative to the disc 
evolves according to the equation 



d No . 2 

— cos p = — sin 
dt Js 



where 



c = 



/?(A-Ccos 2 /3) , (23) 
C'cos 2 6>* 



6rf/2 



(24) 



Equation ( |23[ ) is the key result of this paper. Note that while 
J s depends on the (unspecified) spin-down torque \N S \ , the 
evolution of the spin-orbit inclination angle P depends only 
on two dimensionless parameters: A and £. Although the 
precise values of these two parameters are uncertain (see 
Sect. 2), we expect A to lie in the range of 0.1 - 1, and f to 
range from somewhat less unity to a few (for f ' ~ f ~ 1 and 
?? ~ 0.5). 

Equation ( |23[ ) reveals the following behavior for the evo- 
lution of P (see Fig. [4|): 

(i) For A = 0, equation (23 1 describes the effect of 



the magnetic warping torque acting alone on the star. This 
torque always pushes the stellar spin toward anti-alignment 



5 This expression for f ' assumes that the corotation radius r co is 
larger than r- m t ■ 

6 The coefficient | in equation |2l| is only valid in the limit 
5 = H/r -5- 0. For small but finite 5, we get |[1 - |«/| +0(<5)]. 
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50 100 

/S(degrees) 

Figure 4. The rate of change of the stellar inclination angle j3 
for a fixed disc rotation axis (see Eq, |23| l. From top to bottom, 
the four curves correspond to (A, Q = (1,0.5), (0.8,2), (0.2,2) 
and (0,2), respectively. The arrows indicate the direction of /3 
evolution. For f/A < 1, the stellar spin evolves towards alignment 
for all (see the short-dashed line); for £/A > 1, the spin either 
evolves toward /?+ ^ or toward anti-alignment (/3 = 180°), 
depending on the initial values of /3. The fiducial spin evolution 
timescale t sp i n is defined in Eq. Jl4|. 



with i (see Fig. [I]). In particular, the aligned state (/3 = 0) 
is unstable: For f3 <C 1, we have $//3 — C/^spin- On the 
other hand, the perpendicular state (/3 = 7r/2) represents 
a "bottleneck" where the warping torque J\f w vanishes: For 
P = 7t/2 + A with |A| < 1, we have A/A ~ (C/i spin )A. 
Thus, when the outer disc orientation is fixed (see Section 
5), starting from a small /3, it would take infinite time to 
cross this 90° barrier. 

(ii) For (/X < 1: Regardless of the initial /3, the spin 
always evolves towards alignment. 

(iii) For £/A > 1: There are two possible directions of /3 
evolution, depending on the initial value of /3. The condition 
d cos P/dt = yields two "equilibrium" states (/3+ and /3_), 
given by 



cos^± = ±\/VC 



(25) 



Of the two equilibria, one is stable and the other is 

unstable. For j3{t = 0) < /3_, the system will evolve towards 
a misaligned prograde state for f3(t = 0) > /3_, the 
system will evolve towards the anti-aligned state (/? = 180°). 
Note that the timescale to change the stellar spin 



(Eq. 



14 



can be written as t s 



= (Js/MJ^iM+fM), 



where k n — V GM^r^ is the specific angular momentum 
of the accreting gas at n n . The disk lifetime (observed 
to be around 10 Myrs) is a few times less than M*/M 
if one uses the appropriately averaged value of M. Since 
J s /(M*Zi n ) <C 1, this implies that t ap in is typically less than 
the disk lifetime, i.e., significant change of the stellar spin 
can be achieved during the disk lifetime. Of course, in the 
earlier phases of the protostars (e.g., Class-0 T Tauri, with 
M ~ lO _5 M0/yr), t sp i n is much shorter than 1 Myrs, while 
in the later phases t ap in is longer. 



5 APPLICATION TO SPIN-ORBIT 

MISALIGNMENT IN EXOPLANETARY 
SYSTEMS 

The results of previous sections clearly show that, contrary 
to the standard assumption, the spin axis of a magnetic 
protostar does not necessarily align with the axis of its cir- 
cumstellar disc. Therefore, a planet formed in the disc may 
have its orbital axis misaligned with the stellar spin, even be- 
fore any few-body interactions (such as strong planet-planet 
scatterings and Kozai interactions) take place. A clear pre- 
diction of the expected spin-orbit misalignment angle and 
its distribution for an ensemble of planetary systems is com- 
plicated by the fact that the physics of magnetic star - disc 
interaction is complex: Even with our general parameterized 
model, the evolution of /3 depends on two unknown dimen- 
sionless parameters (A, (). The final misalignment angle and 
its distribution also depend on the usual parameters asso- 
ciated with protoplanetary discs (such as the accretion rate 
and lifetime), as well as on the initial conditions for the disc 
direction relative to the stellar spin. 

We now discuss two scenarios to illustrate the possible 
outcomes for the spin-orbit misalignment angles in exoplan- 
etary systems. 

5.1 Scenario (a) 

Here we consider the case where the collapsing/accreting 
materials that form the protostar and its disc all have an- 
gular momentum axes approximately in the same direction. 
This would lead to an initial state where the stellar spin and 
the disc axis are approximately aligned (with small /3). As 
we showed in Sect. 4, for systems with (/X < 1, the angle /3 
will decrease and the spin-orbit alignment will be enhanced. 
For systems with ("/A > 1, however, the angle f3 will increase 
towards /3+ (see Fig. j4|). 

In the idealized situation where the disc axis I is fixed in 
time, we would expect that some planetary systems to have 
f3 close to zero, while others to have /3 clustered around val- 
ues somewhat less than /3+. Since /3+ = cos -1 \J A/C are 
different for different systems, we may not expect a strong 
clustering. Note that because j3+ is always less than 90°, no 
planet in a retrograde orbit can be produced in this ideal- 
ized situation. However, if we consider the possibility that 
the axis of the outer disc changes in time - as might be ex- 
pected because the gas that feeds the outer disc does not 
necessarily have a fixed rotation axis or because the outer 
disc is perturbed by a distant star in a cluster, then planets 
on retrograde orbits may be produced. To achieve this, the 
outer disc axis must vary with sufficient amplitude [larger 
than (/?_ — and on sufficiently short timescale (shorter 
than the spin evolution time), and the inner disc axis must 
adjust (via bending wave propagation or viscous diffusion) 
quickly to the variation of the outer disc, so that a system 
at the ft < f3+ state may be pushed over to the /3 > /3_ state 
— We will show in Paper II that this is possible. 

5.2 Scenario (b) 

Here we consider the case where the initial stellar spin axis 
and the disc axis are randomly distributed with respect to 
each other. This might be expected if the turbulent gas in 
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Figure 5. Distribution of the angle i/> between the stellar rotation 
axis and the disc axis at t/t apln = (dotted line), 1,2,3,4 (solid 
line), for £ = 0.125 (upper panel) and ( = \f2 (lower panel), both 
with A = 1. The initial distribution is isotropic. 



a molecular cloud core that feeds into the disc falls in with 
random directions. Suppose at the time of planetary for- 
mation, the stellar spin is determined by the accumulative 
angular momentum accretion, then the stellar spin axis may 
be quite different from the "current" disc axis (e.g., Bate et 
al. 2009), and the distribution of the spin-orbit misalignment 
angle may be quite broad. The important problem that we 
wish to address here is: How the distribution function f(ip, t) 
evolves further in time? (Here we define the angle between 
I and U3 3 as ip — /3, in agreement with the notation used in 
observational papers; see Triaud et al. 2010). 

Without the magnetic warping effect discussed in 
Sect. 4, the systems will evolve toward alignment, there- 
fore an initial random distribution, f(ip,0) — sim/j/2, will 
become increasingly more peaked towards ip < 90° as time 
passes. However, when the magnetic warping effect in taken 
into account, a variety of outcomes become possible, as it is 
clear from our results in Sect. 4. 

Consider a simple model where A, £ and t sp i n are con- 
stant in time and the same for different systems. Then the 
evolution of f(ij),t) is governed by the equation 



(26) 



for 



where the time derivative of ip = /3 is given by Eq. ( 23 1 
In Fig. [5] we plot f(ip,t) between t = and t = 4t, 
(A,C) = (1,0.125) and (1,72), assuming an initial isotropic 
distribution f(ip,0) — (sin^)/2. We see that, as expected 
from the behavior depicted in Fig. [4] for £/A < 1, the distri- 
bution function f(ip,t) will become increasingly peaked at 
P < 90°. For C/A > 1, however, f(t[>, t) evolves into a double- 
peaked function, with one of the peaks located around /3+ 
and the other close to 180°. 

Obviously, to determine the true distribution /(i/>, t) for 
an ensemble of planetary systems, one must further "av- 
erage" over the distribution of (A, Q for different systems. 
These parameters are largely unconstrained. Nevertheless, 



we can reasonably expect to see a significant number of sys- 
tems concentrating at small ip (all cases with f/A < 1), and 
the rest distributed between all possible values of /3+ (which 
is less than 90°), or close to ijj = 180°. 



6 DISCUSSION 

We have shown in this paper that the angle between the 
spin axis of a magnetic protostar and the axis of its disc 
may have a wide range of values: some systems are expected 
to be aligned, but some are expected to be misaligned. As 
a result, we would expect the spin-orbit inclination angle 
in exoplanetary systems to have a broad distribution, with 
some systems aligned, while others highly misaligned, even 
without any additional/subsequent physical processes that 
may affect this angle. 

While the results of this paper are compatible with the 
possibility that planet-planet scatterings or Kozai interac- 
tions play a role in determining the orbital characteristics 
of some close-in planets, they nevertheless weaken the chal- 
lenges posed against the disc-migration scenario for the ori- 
gin of most hot Jupiters. Most likely, both few-body interac- 
tions and disc migration are needed to explain the observed 
period distribution of exoplanetary systems. 

Currently, there is no measurement of the spin-disc mis- 
alignment angle in any accreting T Tauri star systems, al- 
though in principle this angle may be constrained from spec- 
tropolarimetry observations (e.g., Hussain et al. 2009; Do- 
nati et al. 2010) and detailed modeling of the variabilities 
of T Tauri stars. Alternatively, spin-disc alignment may be 
tested by measuring the wsinu and the rotation period of 
the star, and comparing the resulting stellar inclination i* 
with the disc inclination idisc, which may be constrained 
from the jet direction (for T Tauri discs) or the sky-projected 
shape of debris disc^] On the theoretical side, the spin-disc 
misalignment angle is determined by the uncertain values 
of two dimensionless parameters, A and C, (see Section 4) 
and their possible dependence on various unknown quanti- 
ties (such as the mass accretion rate) , as well as on the initial 
condition (see Section 5) and the age of the disc-accretion 
phase. However, these two parameters may potentially be 
calibrated from observations of a large sample of T Tauri 
stars. With these input parameters, we can then compare 
theoretical predictions (which include not only the process 
studied in this paper but also few-body interactions that 
may change the spin-orbit angles) with the observed distri- 
bution of misalignment angles between planets' orbits and 
the spins of their mature host stars (Triaud et al. 2010). 

A more direct observational test is to measure the spin- 
orbit angle for multiple planets systems. Few-body inter- 
actions would produce different angles for different plan- 



7 Recently, Watson et al. (2010) carried out such an analysis for 
several debris disc systems and found no significant difference 



between sini* and sinidisc- Note that i* 



does not neces- 



sarily imply alignment between the spin axis and the disc axis. 
Also, systematic uncertainties in estimating idi sc need to be taken 
into account. For example, for HD 22049 (one of the best cases 
studied by Watson et al), the disc inclination is consistent with 
face-on (idi sc fa 25°; K. Stapelfeldt 2010, private communication; 
see Backman et al. 2009). 
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ets in the same system. In contrast, the magnetosphere - 
disc interaction studied in this paper would lead to the re- 
orientation of the stellar spin while preserving the orbital 
plane of the planetary systems. Up to now, there is no 
known multiple-planet systems with measured stellar spin 
and planet orbital angular momentum vector. In the con- 
text of the solar system, the orbits of most planets, including 
those of Jupiter and Saturn, are confined within ~ 2° from 
the ecliptic. Although the Sun's spin vector is misaligned by 
7° with the vector normal to the ecliptic plane of the so- 
lar system, this modest difference may or may not provide 
adequate support for the magnetosphere - disc interaction 
scenario. 

According to the scattering/Kozai scenario, planets 
venture into the proximity of their host stars with a di- 
verse degree of spin-orbit misalignment. Yet, the misaligned 
systems tend to have parent stars with mass greater than 
1.2M (Schlaufman 2010; Winn et al. 2010). Based on 
this scenario, Winn et al. (2010) suggested that spin-orbit 
misalignment for planets around solar type stars has been 
damped out by tidal dissipation in the star, whereas such 
dissipation may be less efficient for more-massive stars (with 
M > 1.2M©) because of their shallow or non-existent outer 
convection zones (e.g., Barker & Ogilvie 2009). Nevertheless, 
the tidal dissipation time scale is a rapidly increasing func- 
tion of planets' orbital semi-major axis. Therefore, according 
to the tidal interaction hypothesis, the spin-orbit misalign- 
ment may nonetheless be preserved for planets around solar 
type stars with periods longer than a few days. 

In the scenario that the spin-orbit misalignment is in- 
duced by the magnetosphere - disc interaction, tidal damp- 
ing of the misalignment angle can still operate more effec- 
tively around solar type stars than more massive stars. In 
addition, the magnitudes of A and £ and therefore the re- 
sultant value of (3 are functions of the stellar mass. In this 
case, the correlation between spin-orbit misalignment and 
the stellar mass would be independent of the planets' orbital 
period (for periods longer than a few days). The confinement 
of multiple planets in a common orbital plane would also 
be preserved. These observational tests, though challenging, 
will eventually provide clues and constraints on the origin of 
hot Jupiters as well as that of the spin - orbit misalignment. 
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